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ABSTRACT 

We analyze the optical spectra of massive (log M^/Mq > 11.4) radio-loud galaxies 
at z ~ 0.2 and z ~ 0.6. Our samples are generated by crossmatching the SDSS 
t^J- \ DR7 and BOSS spectroscopic galaxy catalogues with the FIRST and NVSS radio 

OO ■ continuum surveys. By comparing stellar population parameters of these radio-loud 

I/"") " samples with radio-quiet control samples matched in stellar mass, velocity dispersion 

and redshift, we investigate how the presence of a radio-emitting jet relates to the 
qq ■ recent star formation history of the host galaxy. We also investigate how the emission- 

' line properties of the radio galaxies evolve with redshift by stacking their spectra. Our 

£N) . main results are the following. (1) Both at low and at high redshift, half as many 

radio-loud as radio-quiet galaxies have experienced significant star formation in the 
past Gyr. This difference in star formation history is independent of the luminosity 
of the radio AGN, except at radio luminosities greater than 10 25,5 W Hz -1 , where it 

■ disappears. (2) The Balmer absorption line properties of massive galaxies that have 
^ \ experienced recent star formation show that star formation occurred as a burst in many 

of these systems. (3) Both the radio and the emission-line luminosity of radio AGN 
evolve significantly with redshift. The average [O in] rest equivalent width increases by 
1 dex from z = 0.2 to z = 0.6, and emission line ratios change from LINER-like at low 
redshift to Seyfert-like at high redshift. However, radio galaxies with similar stellar 
population parameters, have similar emission-line properties both at high- and at 
low-redshift. These results suggest that massive galaxies experience cyclical episodes 
of gas accretion, star formation and black hole growth, followed by the production 
of a radio jet that shuts down further activity. The behaviour of galaxies with log 
M*/Mq > 11.4 is the same at z = 0.6 as it is at z — 0.2, except that higher redshift 
galaxies experience more star formation and black hole growth and produce more 
luminous radio jets during each accretion cycle. 
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1 INTRODUCTION 



Large redshi fts surveys such as the Sloan Digital Sky Sur- 
Email: chenym@nju.edu.cn vey (SDSS; lYork et all |200C| 1 and the Two Degree Field 
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(2DF) redshift survey jColless et alj l200ll ) have provided 
the sky coverage needed to define sufficiently large sam- 
ples of nearby radio sources to study the population statis- 
tics and g l obal e nergetics of these systems at 2 ~ 0.1. 
iBest et al.1 |2005ri ) investigated the properties of a sample 
of 2215 nearby radio galaxies with 1.4 GHz radio luminos- 
ity below 10 25 W Hz" 1 from SDSS, finding that the frac- 
tion of radio-loud A GN is a strong function of both black 
hole and stellar mass. iMauch fc Sadler! (|2007l 1 studied a sam- 
ple of 2661 radio-loud AGN selected fr om the 6dF Galaxy 
Survey. They confirmed the findings of I Best et al.1 l|2005bl ) 
that radio-loud AGN preferentially inhabit the brightest and 
most massive host galaxies, and showed that the fraction 
of galaxies which host a radio-loud A GN correlates with 
the infrared isT-band luminosity as L 2 K . IBest et al.1 (|2005bl ) 
also reported that there was no correlation between ra- 
dio luminosity and optical emission-line luminosity for the 
galaxies in their sample, concluding that optical AGN and 
low-luminosity radio-loud AGN are independent phenom- 
ena which are triggered by different physical mechanisms. 
In later work, it was found that at fixed stellar mass, radio- 
loud AGN are preferentially located in denser environments 
than control samples of radio-quiet galaxies with simil ar 
masses and redshifts (IBest et al.ll2007l: Wake et al.1 



Man delbaum et al. 2009 ; [Donoso et al-lboiOll as well as con- 
trol samples of optically-i dentified AGN ( Kauffmann et all 
|200S| ; iDonoso et all 1201(f ). These results have led to the 
paradigm that strong optical AGN activity is associated 
with galaxies with a significant cold gas reservoir and ongo- 
ing star formation (|Heckman et al.ll2004l ). while radio AGN 
activity may be associated with the accretion of hot gas in 

;alaxies located at the centers of mass ive dark matter halos 

Bower et aHl200d ; ICroton et al.ll2006T ). 

Observational support for a scenario in which radio 
AGN play an important role in regulating the growth of 
massive galaxies has accumulated rapidly over the past 
decade. X-ray studies of groups and clusters of galaxies 
with Chandra and XMM-Newton have shown that these 
jets interact strongly with their environment, blowing clear 
cavities or 'bubbles' in the surrounding X-ray-emitting gas 
jBoehringer et al]|l993l: iFabian et al.ll2000l; IChurazov et al.l 
120011 ; IFabian et al.1 12003 . 120051 ; iForman et al.ll2005l ). These 
cavities provide a d irect estimate of the radio jet power 
jBirzan et al.1 |2004 | 20081 ) and hence an empirical cal- 
ibration of the relation between radio luminosity and 
mechanical power of the radio jet in nearby galaxies. 
The energy input from jets has been found to be suffi- 
cient to prevent star formation activity in massive early- 
type galaxies by heating the int erstellar gas and sup- 
pressing the onset of cooling flows jBinnev &: Tabor] Il995| ; 
Birza n et aUl2004lRawlings fc Jarvisll2004l ; IBest et aUl2006 



20071 ; ISchawinski et al.ll2006l ). 

If radio jets regulate ongoing star formation in massive 
galaxies, one might expect to see differences in the stellar 
populations of the host gala xies of radio-loud AGN com- 
pared to radio-quiet galaxies. IBest et al.1 ((2005b) compared 
the fraction of radio-loud AGN among galaxies with differ- 
ent values of the 4000 A break strength (D4000) , which is a 
powerful measure of the mean age of the stellar population in 
a galaxy. As can be seen from Figure 10 of that paper, there 
appears to be a lower fraction of radio-loud AGN among 
the most massive galaxies (> 1O 11 M0) with young stellar 



populations (i.e. low values of D4000). The statistics were 
rather poor, so Best et al. did not place much emphasis on 
this result. 



While much of our knowledge about radio AGN is based 
upon studies of local galaxies, our understanding of the 
interplay between gas cooling processes and radio AGN 
feedback at higher redshifts is still very much in its in- 
fancy. Strong cosmological evolution of the high-luminosity 
radio source population has been long established (e.g., 
iDunlop "fc^ Pcacock 1990), but the role of this strongly evolv- 
ing population in regulating the growth of massive galaxies 
at hig h redshifts is largely unconstrained. Ijohnston et al.l 
(2008) constructed high signal-to-noise ratio (S/N) com- 
posite optical spectra of 375 radio-detected luminous red 
galaxies at 0.4 < 2 < 0.8. They compared composite spec- 
tra binned by radio luminosity and redshift with those of 
radio-quiet control galaxies and concluded that there were 
no differences in the stellar populations of the two sam- 
ples, except for the most powerful radio galaxies (with 
ii.4GHz > 10 26 W Hz -1 ), which had stronger optical emis- 
sion lines and younger stellar populations. We note, however, 
that the galaxy survey on which this analysis was based was 
subject to a strong pre-selection on galaxy colour and would 
miss blue, star-forming galaxies enitirely. 



The Baryon Oscillation Spectroscopic Survey (BOSS; 
Da wson et al. 2012, submit ted), one of the surveys of SDSS- 
III (|Eisenstein et al]|201ll ). is currently obtaining spectra 
of 1.5 million luminous galaxies to 2 ~ 0.7. At redshifts 
greater than 2 ~ 0.55, the galaxy samples are reasonably 
complete at stellar masses greater than 10 n ' 4 Mo. In this 
paper, we have cross-correlated galaxies with 2 > 0.55 from 
the BOSS spectroscopic survey with the National Radio As- 
tronomy Observ atory (NRAO) Very Large Array (VLA) Sky 
Survey fNVSS; ICondon et al.lll998l ) and the Faint Images 
of the Radio Sky at Twenty centimeters (FIRST) survey 
l|Becker et al.|[l995T ) and have generated the largest sample 
of intermediate redshift radio galaxies to date. We use this 
sample to investigate whether radio-loud AGN have had un- 
usual star formation histories (SFHs) com pared to galax- 
ies of similar mass and redshift. Similar to Ijohnston et all 
(2008), we create composite optical spectra to study the 
emission line properties of radio-loud and radio-quiet galax- 
ies. By comparing results obtained for the BOSS samples 
with those from the lo w-redshift radio AGN catalogue of 
IBest fc Heckmanl (|2012l ). we study how the radio AGN pop- 
ulation has evolved over a timescale of 3.3 Gyr. 



Our paper is arranged as follows. In §2, we describe the 
low and high redshift radio-loud AGN samples and the con- 
struction of the radio-quiet control samples. In §3, we dis- 
cuss how we extract information about the past star forma- 
tion histories of the galaxies in our samples, present results 
for the radio-loud AGN hosts and the control galaxies, and 
analyze the emission line properties of radio galaxies using 
composite spectra . Our results are summarized in §4. We 
use the cosmological parameters Ho = 70 km s _1 Mpc" 1 , 
SIm = 0.3 and Qa = 0.7 throughout this paper. 
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Figure 1. The distribution of the stellar mass and radio luminosity for the DR7 (red) and CMASS (black) radio- loud AGN samples. 
Our DR7 sample includes 5818 galaxies with 0.1 < z < 0.3, while the CMASS sample includes 2856 galaxies with 0.55 < z < 0.7. 



2 SAMPLE SELECTION 

2.1 The low-redshift radio-loud AGN sample 

The Sloan Digital Sky Survey (SDSS; lYork et alJl200Ch ob- 
tained photometry of nearly a quarter of the sk y. A drift- 
scanning mosaic CCD camera (|Gunn et alJll998h mounted 
on the SPSS 2.5m telescope at Apache Point Observatory 
|Gunn et al.ll2006t) imaged the sky in the u, g, r, i, z bands 
(Fukugita ct al. 1996). The imaging data are astrometrically 
dPier et all l2003f) and ph otometrically (|Hogg et alj l200ll ; 
Padmana bhan et al. 2008) calibrated and used to select 
stars, galaxies, and quasars for follo w-up fiber spectroscopy . 

The seventh data release (DR7; lAbazaiian et alj 12009ft 
of the SDSS includes about one million galaxy spectra. The 
spectra have a wavelength coverage of 3800 - 9200A and 
are taken through 3" diameter fibers. The low-redshift ra- 
dio galaxy sample is described in iBest fc Heckmanl (|2012l ). It 
has a flux density limit of 5mJy, which means the sample ex- 
tends down to a radio luminosity of Li.4GHz « 10 23 W Hz" 1 
at redshift z — 0.1. This sample was constructed by 
cross-correlating 927,552 galaxies from the DR7 MPA/JHU 
value-added spectroscopic catalogue^] with the NVSS and 
FIRST surv eys. The cr o ss-mat ching followed the method 
described by IBest et all (l2005al). but also in cluded the im- 
provement developed bv lDonoso et al.l (l2009h for identifica- 
tion o f sources without FIRST detections. Best fc Heckmanl 
(2012) also separated radio AGN from star-forming galaxies. 
We refer the reader to their Appendix A for more details. 

For the current paper, we begin with a subset of 5818 
radio galaxies with r-band magnitude 14.5 < r < 17.77, 
reds hift 0.1 < z < 0.3 and SP ECPRIMARY flag = 1 drawn 
from lBest fc Heckmanl (|2012l ) catalogue. The r-band magni- 
tude limit restricts t he sources to the 'main' galaxy sample 
l|Strauss et alj|2002h . With the SPECPRIMARY flag = 1, 
we avoid repeat observations of a given object. Properties of 
the radio source host galaxies required in this work include 



1 The MPA/JHU catalogue i s 
http://www.mpa-garching.mpg.de/SDSS/ 



available at 



total stellar masses (M*), stellar velocity dispersions (ct») 
and SFH parameters de rived from the pr incipal component 
analysis method (PCA. IChen et al.ll2012T ). Radio luminosi- 
ties are calculated using the formula 



log 10 [Li.4GHz] = log 10 [47rZ)£(z)Sl.4GHz(l + z) 



(1) 



where Dl is the luminosity distance in the adopted cosmol- 
ogy, 5i.4GHz is the measured radio flux density, (l + z)"^ 1 is 
the standard fe-correction used in radio astronomy and a is 
the radio spectral index (S v ~ v~ a )- We ad opted a = 0.7, 
as is usually assumed for FRI radio galaxies (jCondon et al.l 
2002). The red histograms in Figure \T\ show the distribution 
of stellar mass and radio luminosities of the low redshift 
radio-loud AGN sample. Any galaxy that meets our radio 
flux cut is termed "radio-loud" in this paper. 

2.2 The high- redshift radio-loud AGN sample 

The SDSS-III project has completed an additional 3000 deg 2 
of imaging in the southern Galacti c cap in a manner i denti- 
cal to the original SDSS imaging (jAihara et alj|201ll ). The 
Data Release 9 (DR9) of BOSS includes spectra of 462,979 
luminous galaxies to z ~ 0.7 (Ahn et al. 2012, submitted). 
The spectrographs (Smee et al. 2012, submitted) have been 
significantly upgraded from those used by SDSS-I/II, with 
improved CCDs with better red response, high throughput 
gratings. The new fibers are 2" in diameter and the spectra 
cover the wavelength range 3600-10,000A. Galaxy spectro- 
scopic redshift is determined by the BOSS pipeline (Bolton 
et al. 2012, submitted). The details of ho w targets are se- 
lected from the photometry are described in lEisenstein et ahl 

The high redshift galaxy sample we analyze here is the 
"CMASS" sample (so-named, because it is very approxi- 
mately stellar-mass limited). This sample is defined by the 
following cuts: 



d± > 0.55 and 17.5 < i < 19.9 and i flbro 2 < 21.5 



(2) 



i < 19.86 + 1.6 (dj. - 0.8) and r - 1 < 2 
where d± is a "rotated" combination of colours defined as 
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d± — (r — i) — (g — r)/8. The quantity Zfibrc2 is the i-band 
magnitude of the galaxy measured wit hin a 2" BOSS fibre i n 
the SDSS ugriz photometric system dFukugita et al.lll996l ). 
All magnitude limits are given in terms of "cmodel" magni- 
tudes whereas colour cuts are defined using "model" mag- 
nitudes. Two additional conditions are introduced to reduce 
contamination by stars: z ps { — z mo doi ^ 9.125 — 0.46z mo d e i 
and i pS f - i mo dci > 0.2 + 0.2 x (20.0 - i mo dd) ("psf" refers 
to the psf Mag quantity in the SDSS database). 

The cuts listed in equation (2) are designed to identify 
massive galaxies at z > 0.4. The parameter d± > 0.55 se- 
lects galaxies at high enough redshift that the 4000 A break 
has shifted beyond the observer frame r-band, leading to 
red observed r — i colours. The cut on the i-band magni- 
tude is designed to produce a sample that is approximately 
complete down to a limiting stellar mass of M* ~ 10 n ' 2 MQ 
(|Chen et alj|2012h . 

We cross match the CMASS galaxie s with NVSS and 
FIRST , using the same method as in iBest fc Heckmanl 
(2012). We generate a sample of 5469 radio galaxies with a 
flux density limit of 3mJy. In this work, we limit the sample 
to 2856 galaxies with z > 0.55 where the colour incomplete- 
ness of the sampl e is minimized. As discussed in detail in 
I Chen et~aL I <|2012l ). the > 0.55 colour cut used to define 
the CMASS galaxies results in blue galaxies b eing missed 
from t he survey, particularly at lower redshifts. IChen et al.l 
(|2012h showed that the fraction of massive galaxies with re- 
cent star formation could be under-estimated by more than 
an order of magnitude at z < 0.55. At higher redshifts, the 
correction to the fraction of actively star-forming galaxies is 
approximately a factor of two. The 3mJy flux limit corre- 
sponds to a radio luminosity of Li.4GHz ~ 10 24,5 W Hz -1 at 
redshift z = 0.55. At these very high radio luminosities, we 
will almost exclusivel y pick up radio-loud AGN rather tha n 
star-forming galaxies (| Sadler et alj|2002l ; [Best et al.ll2005bh . 

Once again, the properties of the CMASS radio source 
host galaxies, including total stellar masses, stellar ve- 
locity dispersions and SFH parameters, are derived from 
the principa l com ponent analysis method, as described in 
l|Chen et al.l I2012T ). These authors found that the PC A 
method gives more robust estimations of D4000, H<5A and 
a, for the low S/N CMASS spectra (> 85% galaxies have a 
median S/N per pixel (69km s -1 ) below 4). The black his- 
tograms in Figure[T]show the distribution of stellar mass and 
radio luminosity of the high redshift radio-loud AGN sam- 
ple. As can be seen, the high-redshift radio AGN are offset 
to higher stellar masses and radio luminosities compared to 
the low-redshift radio AGN, but there is a reasonable range 
in mass and luminosity where direct comparisons between 
the two samples can be made. 



2.3 Radio-quiet control samples 

We create samples of radio-quiet galaxies, which are 
matched as closely as possible to the radio-loud AGN hosts. 
The purpose of the control samples is to determine which 
galaxy properties are causally linked with the radio-emitting 
jet. 

For each galaxy in the radio-loud sample, we find an- 
other galaxy located in the FIRST and NVSS survey area 
without a radio detection that is closely matched in redshift 



(|Ajz| < 0.005), stellar mass (|AlogM*| < 0.05), and velocity 
dispersion (|A<r.| < 10km s" 1 ). 

The motivation for choosing this set of matching pa- 
rameters is the following: 

(i) By constraining the control galaxy to be at similar 
redshift, we avoid any confusion due to evolutionary effects 
and we make sure the SDSS fibre spectrum includes light 
coming from the the same physical scales in both the radio 
galaxy and the control. 

(ii) Constraining the control galaxies to have similar stel- 
lar mass is extremely important because the probability that 
a galaxy is a radio source is a strong function of mass. Stel- 
lar population properties are known to vary strongly as a 
function of galaxy mass. 

(iii) Constraining the control galaxies to have similar ve- 
locity dispersions ensures that the radio-loud and the radio- 
quiet galaxy samples have roughly the same black hole mass 
distributions. Any AGN physics that is determined primar- 
ily by the black hole mass will thus be the same for the two 
samples. 

Any radio galaxy without a radio-quiet twin is excluded from 
our analysis. This leaves a total of 5416 (2820) radio galaxies 
in the DR7 (CMASS) samples. 



3 RESULTS 

3.1 Stellar populations of radio-loud and 
radio-quiet galaxies 

In this section, we analyze the standard absorption line in- 
dices D4000 and H<5A, reconstructed using our PCA tech- 
nique, as well as the parameter F t , defined as the ratio be- 
tween the mass of the stars formed in the last Gyr and the 
stellar mass of the galaxy. The quantity F* can be trans- 
formed into the star formation rate (SFR) averaged over the 
last Gyr using the equation, SFR = F, x A'/«/10 9 M yr -1 . 

Figure [2] presents histograms of these quantities for 
DR7 radio-loud galaxies and controls in the top panel, and 
for CMASS radio-loud galaxies and controls in the bottom 
panel. The radio-loud galaxy sample (black histograms) is 
clearly shifted to larger values of D4000, smaller values of 
H<5A and smaller values of The largest differences are 
seen for H5A and for F, at both low and at high red- 
shifts. We performed a Kolmogorov-Smirnov (K-S) test on 
the D4000, H<5A, F t distributions of radio-loud and radio- 
quiet samples and obtained a probability smaller than 10 -5 
for these three parameters for the null assumption that the 
radio-loud and radio-quiet galaxies are drawn from the same 
parent distribution. 

We now focus on (1) quantifying the difference between 
the stellar populations of radio-loud and radio-quiet galaxies 
in more detail; (2) studying the evolution of this difference 
with redshift. We define a parameter R = Nq/Nl, where 
Nl (Nq) is the number of galaxies with F, > F* iCr it in the 
radio-loud (quiet) samples. F* tCr i t is a certain value of F*. 
The left panel of Figure [3] shows J? as a function of logF* jCr i t 

2 The number of galaxies included in the F* histograms is smaller 
than in the other two, because the majority of galaxies have un- 
detectable small F„. 
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Figure 2. Comparison of galaxy stellar population parameters D4000, H<5A and F» for the radio-loud galaxy samples and their radio- 
quiet control samples. The top panel shows results for DR7; the bottom panel for CMASS. The radio-loud galaxies are shown as black 
histograms while the radio-quiet galaxies are shown as red histograms. The radio-loud and radio-quiet samples include exactly the same 
number of galaxies. There are 5416 and 2820 galaxies in the DR7 and CMASS samples, respectively. 




Figure 3. The ratio R between the number of galaxies with F, > F* c rU m the radio-loud and radio-quiet samples is plotted as a 
function of logF* cr £j. In the left panel, we show results for DR7, while right panel is for CMASS. In each panel, the red line represents 
galaxies with Fi.4GHz larger than the median radio luminosity (10 24 W Hz -1 for DR7, 10 25 W Hz -1 for CMASS), while the blue line is 
for galaxies with Fi.4GHz smaller than these values. The black line is the result for the whole sample. The diamonds in the bottom-left 
corner indicate the values of R for galaxies with undetectable small F* . 



for DR7 galaxies, while the right panel shows the same thing 
for CMASS galaxies. In each panel, the red line represents 
galaxies with Liaghz larger than the median radio luminos- 
ity (10 24 W Hz" 1 for DR7, 10 25 W Hz" 1 for CMASS), while 
the blue line is for galaxies with Liaghz smaller than these 



values. The black line represents the whole sample. The di- 
amonds in the bottom-left corner represent the ratio be- 
tween the numbers of radio-quiet and radio-loud galaxies 
with F, ~ 0. Figure indicates that, 1) the ratio R is al- 
ways larger than 1, implying that the ratio of star forming 
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Figure 4. We match the DR7 and CMASS radio-loud samples by stellar masses. Left panel: the distributions of logF* for the mass- 
matched samples, solid lines for DR7, dashed lines for CMASS, black and red lines colour-code the radio-loud and radio-quiet galaxies, 
respectively. In both DR7 and CMASS, there are more actively star-forming galaxies with logF, > —2 (the vertical blue line) in the 
radio-quiet sample than in the radio- loud sample. Similar to Figure [3] the right panel shows R as a function of -F tjCr ;t for the DR7 
(black) and CMASS (red) samples. 



to passive galaxies is always higher in the radio-quiet sam- 
ple, 2) the ratio R increases with F,. cr i t , implying that radio 
AGN are more important for suppressing strong star forma- 
tion as compared to weak star formation, 3) the differences 
in R for low-luminosity and high-luminosity radio galaxies 
are negligible. 

From Figure [31 one might be tempted to infer that 7? is 
larger at higher redshifts. However, it is not fair to compare 
the DR7 and CMASS samples, because these two samples 
have very different stellar mass distributions (see the left 
panel of Figure[T]), and the star formation history of a galaxy 
strongly depends on its stellar mass (e.g. , IZheng et alTl2007l : 
IChen et al]|2009l ; iKarim et al.ll201ll ). In order to study evo- 
lutionary effects, we match the DR7 and CMASS galaxies 
by stellar mass. We bin both samples by stellar mass with 
|AlogM*| < 0.02. For each stellar mass bin, we count the 
number of galaxies in both DR7 (nl) and CMASS (n2) sam- 
ples. If nl > n2, we randomly select n2 galaxies from this 
certain mass bin of DR7, otherwise, we select nl galaxies 
from CMASS. This matching process leaves 1658 galaxies 
in both DR7 and CMASS samples. Unless specified other- 
wise, all future plots are based on mass-matched samples. 

The left panel of Figure U shows the distributions of 
logF* for the four mass-matched samples: solid lines show 
results for DR7, dashed lines for CMASS; distributions for 
radio-loud and radio-quiet galaxies are shown as black and 
red histograms, respectively. Two conclusions are evident 
from this panel: (1) the distribution of F* is shifted to much 
higher values in the CMASS sam ple than in the DR7 sample. 
As shown bv lChen et aI1l|2012l ). star formation in very mas- 
sive galaxies (M„ > 1O 11,4 M0) has evolved strongly since a 
redshift of ~0.6; (2) in both DR7 and CMASS, there are 
more actively star-forming galaxies with logF* > —2 (the 
vertical blue line) in the radio-quiet sample than in the 
radio-loud sample. The right panel shows -R as a function 
of F* iCr it for the mass-matched DR7 (black) and CMASS 
(red) samples. As can be seen, the range in R spanned by 



the two samples is much the same. For both the DR7 and 
CMASS samples R increases from ~ 1 at low F, tCr it to a 
maximum value of around 2 — 2.5. The main difference is 
that the distribution of F* has been shifted to higher values 
in the CMASS sample, so the maximum occurs at a higher 
value of F ttCr it. 

We now ask whether the fraction of actively star form- 
ing galaxies in the radio-loud sample depends on the thresh- 
old radio luminosity used to define it. We set logF„ iCT -i t = 
-1.5 (-2) for the CMASS (DR7) samples and we define the 
parameter <j> = N^t/N^n to represent the fraction of galax- 
ies with f* > F, tCr i t . Figure [S] shows <f> as a function radio 
luminosity. Black lines are for DR7, while red lines are for 
CMASS. The dashed and solid lines represent radio-loud and 
radio-quiet galaxies, respectively. As expected, the fraction 
of actively star forming galaxies is higher in the radio-quiet 
samples. The main point of this plot is that <f> remains con- 
stant over a wide range in radio luminosity. At the very high- 
est radio luminosities (Li. 4 ghz > 10 25 5 W Hz -1 for DR7, 
il.4GHz > 10 26 W Hz" 1 for CMASS), <j> suddenly increases 
and the difference in the fraction of actively star forming 
galaxies in the radio-loud and the radio-quiet samples dis- 
appears. 



3.2 Constraints on the star formation histories of 
massive galaxies 

Figure [2] shows that the difference in H<5A between radio- 
loud AGN and radio-quiet control galaxies is much larger 
than the difference in D4000. This result i s also true for the 
mass- m atched samples. As discussed in Ka uffmann et af] 
(|2003al lbl), the location of galaxies in the D4000— H<5A plane 
can tell us whether they have been forming stars in a burst or 
continuously over the last 1 — 2 Gyr. Galaxies with continu- 
ous star formation occupy a narrow strip in the D4000— H<5A 
plane, while a recent burst displaces them away from this 



© 0000 RAS, MNRAS 000, 000-000 



The link between radio AGN activity and star formation 7 




Figure 6. The DR7 (left) and CMASS (right) galaxies on the D4000— H<5A plane. Black dots represent radio-loud galaxies and grey dots 
represent radio-quiet galaxies. The over-plotted green line shows the locus of a model galaxy with exponentially declining SFH with time 
constant r = 2 Gyr and twice solar metallicity. The orange-cyan line shows how a galaxy with an underlying old stellar population that 
experiences a recent burst of star formation evolves in the D4000— H<5A plane. This model is calculated using Bruzual & Chariot (2003) 
population synthesis templates. The SFH consists of two parts: an underlying continuous model with r = 1 Gyr plus a recent burst. The 
burst, which has twice solar metallicity and duration 0.5 Gyr, occurs when the galaxy is 8 (5) Gyr old and the fraction of stellar mass 
produced during the burst relative to the total mass formed by the continuous model is ~1% (~3%) in the left (right) panels. The orange 
diamonds mark the place where the burst begins. Following the burst, every colour-coded section of the track represents a duration of 
100 Myr. The median ±1(7 measurement errors on D4000 and H<5A are indicated in the right side of each panel. 




23.5 24.0 24.5 25.0 25.5 26.0 26.5 

lOg L 1.4GHz W Hz" 1 

Figure 5. This figure shows the fraction of actively star form- 
ing galaxies as a function radio luminosity. Actively star-forming 
galaxies are defined to be those with logF» iCrit > —1.5 (—2) for 
the CMASS (DR7) samples. Black lines show results for DR7, 
while red lines are for CMASS. The dashed and solid lines repre- 
sent the radio-loud galaxies and the radio-quiet controls, respec- 
tively. Note that for the radio-quiet galaxies, we use the radio 
luminosity from their radio-loud twins as the x-axis quantity. RL 
and RQ represent radio-loud and radio-quiet respectively. 



locus to higher values of H5A for ~ 1 — 2 Gyr following the 
burst. 

Figure E] shows the DR7 (left) and CMASS (right) 
galaxies in the D4000— H5A plane. Black dots represent 
radio-loud AGN and grey dots represent radio-quiet galax- 
ies. The over-plotted green line shows the locus of a model 
galaxy with exponentially declining SFH with time constant 



t = 2 Gyr and twice solar metallicity. The orange-cyan line 
shows how a galaxy with an underlying old stellar pop- 
ulation that experiences a recent burst of star formation 
evolv es in the D4000— Ht5A plan e. This model is calculated 
using iBruzual fc Chariot! (|2003l ) population synthesis tem- 
plates. The SFH consists of two parts: an underlying contin- 
uous model with r = 1 Gyr plus a recent burst. The burst, 
which has twice solar metallicity and duration 0.5 Gyr, oc- 
curs when the galaxy is 8 (5) Gyr old and the fraction of 
stellar mass produced during the burst relative to the total 
mass formed by the continuous model is ~1% (~3%) in the 
left (right) panels. In the plot, the orange diamonds mark 
the place where the burst begins. Following the burst, every 
colour-coded section of the track represents a duration of 
100 Myr. 

These models are meant to be purely illustrative, and 
the values of burst duration and mass fraction that we have 
adopted are not tuned to fit to the data. The main con- 
clusion that one reaches from studying Figure [6] is that the 
majority of CMASS and DR7 galaxies with D4000 < 1.7 
lie systematically above the locus occupied by galaxies that 
have experienced exponentially declining star formation his- 
tories. These galaxies have D4000/H5A values that can be 
explained if they have experienced recent bursts of star for- 
mation. 

There are many more galaxies in the post-burst phase in 
which H<5A is offset to large values than in the "young burst 
phase", when H5A is low (in fact, often in emission). Part 
of the reason for this behaviour is that the time spent in the 
young burst phase is short. In addition, we caution that the 
d± colour cut used to define CMASS sample may exclude 
galaxies with a significant population of very young stars. 
Figure[7]shows the track of the burst model in the right panel 
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Figure 7. The predicted behaviour of the colours of the burst 
model in the right panel of Figure l6l with M* = W 11 5 Mq at 
z = 0.6. This figure shows the locus of this burst model in the 
dj_ — i plane, where i represents the i-band magnitude. The or- 
ange diamond marks the point when burst begins. Following the 
burst, every colour-coded part of the locus represents a duration 
of lOOMyr. The horizontal dashed line marks the CMASS d ± cut. 



of Figure [6] in the dj_—i plane, where i is the i-band magni- 
tude of the galaxy. A galaxy with a stellar mass of 10 Mq 
at the end of the burst is assumed. Once again, the orange 
diamond marks the point when burst begins; following the 
burst, every colour-coded part of the locus represents a du- 
ration of lOOMyr. The horizontal dashed line marks the d± 
cut. It is clear that a 3% burst can result in blue enough 
colours that a fraction of galaxies at z = 0.6 may drop out 
of the sample. The burst track plotted in Figure [7] assumes 
no observational errors, which will scatter some young bursts 
into the region of d± > 0.55. The effect of dust on the colours 
has also not been included in the models. 

3.3 Emission line properties from composite 
spectra 

3. 3. 1 Stacking by stellar population age 

In this section, we analyze composite spectra in bins of F*, 
D4000 and USA for the DR7 and CMASS samples. The 
galaxy spectra are first corrected fo r foreground Galacti c 
attenuation using the dust maps of ISchlegel et alj (|l998h , 
transformed from flux densities to luminosity densities, and 
shifted to the rest frame using the redshift determined by 
the SDSS DR7/BOSS pipeline. Since most of the galaxies 
in our samples are intrinsically luminous galaxies, we nor- 
malize each spectrum by its mean luminosity in the wave- 
length range 5300-5500A, and simply average the spectra 
with equal weight (note that the weight of the bad pixels 
identified in the SDSS mask array is set to zero). 

Figure[8]shows normalized composites of CMASS galax- 
ies with logF, < -2 (black), -2 < logF, < -1.5 (red) and 
logF* > —1.5 (blue). We have not attempted to divide the 
galaxies with logF, < —2 into separate bins because the 
PCA method is not able to recover the SFHs of galaxies in 
which less than ~1% of the stellar population formed in the 
last Gyr. As expected, for both radio-loud (left panel) and 



radio-quiet (right panel) galaxies, emission line fluxes (e.g. 
[O n], [O ill]) increase for larger values of F* and the spec- 
tra also become bluer, and the bluest galaxies have stronger 
emission lines if they are radio- loud than radio-quiet. Com- 
posites binned by D4000 and H5A yield very similar results, 
so we do not present them here. 



3.3.2 Stacking by radio luminosity 

Figure [9] shows the difference between the composite spec- 
tra of radio-loud and radio-quiet galaxies for different bins 
in radio luminosity. Radio luminosity increases from top to 
bottom in the plot. The left panel shows results for DR7 and 
the right panel for CMASS. The red lines mark zero flux 
density difference. The near zero residual flux in Figure |9j 
except for regions around the emission lines, indicates that 
the average stellar continua of radio-loud and radio-quiet 
galaxies are quite similar. At first sight, this result may ap- 
pear to conflict with the conclusions pre sented in the previ- 
ous se ction. However, when we check the lBruzual fc Charlotl 
(2003) stellar population models, we find that a simple stel- 
lar population with an age of 0.5 Gyr is typically 10 times 
less luminous in the r-band than a 5 Gyr old stellar pop- 
ulation of similar mass. In the CMASS radio-quiet sample, 
the fraction of galaxies that have formed more than ~3% of 
their stellar mass in the last Gyr is about 10%. We estimate 
that if we stack all the CMASS radio-quiet galaxies together, 
the fraction of light contributed by stars less than 1 Gyr in 
age will contribute only ~ 3% of the total flux in the com- 
posite spectrum. This explains why the stellar continua of 
the radio-loud and radio-quiet composite spectra are nearly 
identical. 

There are almost no emission lines visible in the 
composite spectra of the DR7 galaxies, except weak 
[O in] AA4959,5007, which appears in panel (d) in the com- 
posite spectra of the most powerful radio sources. Almost 
all very massive galaxies (M* > 10 11a Mq) are inactive at 
low redshifts. In contrast, the spectra of the CMASS radio- 
loud galaxies have obvious emission lines (e.g., [O n], [Ne ill], 
H/3, [O ill]) and the emission line strengths increase strongly 
with radio power. 

There are two bins of radio luminosity where we can 
compare the DR7 and CMASS composite spectra directly. 
Figure [10] shows the residual between CMASS and DR7 
composites for these two bins: black lines are for radio- 
loud galaxies, while red lines are for radio-quiet galaxies. 
It is clear that the CMASS composite spectra are bluer and 
have stronger Balmer absorption lines. This indicates that 
the high redshift massive galaxies have significantly younger 
stellar populations on average. 

3.3.3 Correlations between emission line properties and 
star formation 

Follow ing the methodology describ ed in iTremonti et al.l 
l|2004l) and lBrinchmann et al.l l|2004T ). we fit a linear super- 
position of stellar population templates to the stellar con- 
tinua of the composite spectra presented in §3.3.1 & 3.3.2. 
We subtract the continuum model from the composite, and 
measure the rest equivalent width (EQW) of [O iii]A5007 
from the residual spectrum. 
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Rest-frame Wavelength (A) 

Figure 8. Composite spectra of CMASS galaxies in different F* bins. The left panel shows results for radio-loud galaxies and the right 
panel for radio-quiet galaxies. RL and RQ represent radio-loud and radio-quiet respectively. 
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Figure 9. The difference between the composite spectra of radio-loud and radio-quiet galaxies for the different radio luminosity bins. 
Radio luminosity increases from top to bottom. The left panel shows results for DR7 and the right panel for CMASS. The red lines mark 
the zero point (i.e. no difference between the spectra). 
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Figure 10. The difference between CMASS (high-z) and DR7 (low-z) composites. Black lines show results for radio-loud galaxies, while 
red lines show results for radio-quiet galaxies. It is clear that the composites of CMASS galaxies are bluer and have stronger Blamer 
absorption lines than the DR7 galaxies although they have similar stellar masses. 
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Figure 11. This figure shows the correlations between [O in] EQW and D4000 (left), radio power (right) for the composites in various bins 
of D4000 and radio power. Black symbols represent radio-loud galaxies and red symbols represent radio-quiet galaxies. Stars connected 
by dashed lines represent galaxies from DR7 and solid points represent galaxies from CMASS. We estimate the sample variance errors 
using a bootstrap technique. RL and RQ represent radio-loud and radio-quiet respectively. 



Correlations between [O ill] EQW and D4000 and ra- 
dio power are shown in Figure [11] Black and red sym- 
bols represent radio-loud galaxies and radio-quiet galax- 
ies, respectively. Stars connected by dashed lines represent 
galaxies from DR7 and solid points represent galaxies from 
CMASS. As can be seen from the left panel, the strength 
of [O in] emission correlates strongly with D4000 for both 
radio-quiet galaxies and for radio-loud AGN. In the DR7 
galaxy sample, the emission lines from extremely massive 
galaxies (> 10 114 Mq) almost always have ratios indicative 
of ionization by "hard spectr um" sources (AGN or post- 
asymptotic giant branch stars; Kauffma nn et al. | |2003<J ). For 
the CMASS galaxies, [N n] and Ha have shifted out of 
the observed spectral range, but alternative diagnostic di- 
agrams such as E QW-[Q n]/EQW-H/3 vs. log EQW-[0 11] 
|Rolaet al.lll997l ) indicate that the emission is still domi- 
nated by AGN. 

If we regard the [O in] line as a tracer of the accre- 



tion of matter onto a central black holfl the left panel of 
Figure[TT]indicates that the average black hole growth is cor- 
related wit h the growth of the ga laxy, as has been discussed 
previously |Heckman et al.ll2004h . It is interesting that at 
fixed D4000, radio AGN have somewhat higher [O in] EQW 
values than the control galaxies. This suggests that shock 
excitation from the jet may contribute to the [O ill] emis- 
sion. The effect is quite weak (less than a 50% change in 
[O in] EQW). To zeroth order, one relation between [O ill] 
EQW and D4000 seems to adequately describe radio-loud 
and radio-quiet galaxies at z ~ 0.2 and at z ~ 0.6. 

The right panel of Figure [TT] shows trends in [O ill] 
EQW as a function of radio power for the composite spec- 
tra constructed from the DR7 and CMASS samples. This 

3 Population synthesis models (e.g. Bruzual &: Chariot 2003) 
show that [OIII] EQW of more than a few A cannot be produced 
by PAGB star ionization. 
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Figure 12. The line ratio between [O III] and H/3 as a function of D4000 (left panel) and radio power (right panel), the symbols and 
lines have the same meaning as that in Figure II II The blue ho rizontal line demarcates the value of [O III] /H/3 that separates Seyferts 
from LINERs in the AGN population llKauffmann et al.|[20Q3c]) . We estimate the sample variance errors using a bootstrap technique. 
RL and RQ represent radio-loud and radio-quiet respectively. 



plot shows that the [O ill] EQW does not correlate with ra- 
dio power at radio luminosities below ~ 10 25 ' 5 W Hz" 1 both 
for the DR7 and the CMASS radio samples. At fixed radio 
luminosity, high-redshift radio galaxies have higher [O ill] 
EQW than low-redshift radio galaxies. The same effect is 
seen in the control sample, so we attribute this evolution to 
the fact the amount of recent star formation in the massive 
galaxy population increases with redshift. At radio luminosi- 
ties higher than ~ 10 25 ' 5 W Hz" 1 , the [O ill] EQW increases 
as a function of radio power. This is not seen in the control 
sample galaxies, so this phenomenon must be intrinsic to the 
radio-emitting jet itself. 

Figure [12] shows the ratio between the [O ill] and H/3 
line fluxes as a function of D4000 and radio power. The 
symbols and lines have the same meaning as that in Fig- 
ure 1111 The blue horizontal line demarcates the value of 
[O in] /H/3 t hat separates Seyferts fro m LINERs in the AGN 
population l|Kauffmann et all 12003d ) . The trends visible in 
this plot mirror what is seen in Figure [TT1 There is a strong 
correlation between [O m]/H/3 and D4000, which does not 
depend on redshift and which is roughly the same for radio- 
loud and radio-quiet galaxies (we again see a slight excess in 
[O ill] emission for the radio-loud population in this plot; it 
is particularly evident at low values of D4000). The "univer- 
sal" correlation between [O in] /H/3 and D4000 again leads 
us to attribute much of the evolution in [O ill] /H/3 in the 
right panel of Figure [12] to the changing stellar populations 
of high-redshift galaxies. However, the abrupt increase in 
[O ill] /H/3 to Seyfert-like values found at radio luminosities 
m excess of ~ 10 25 ' 5 W Hz" 1 is not seen in control samples, 
and must be intrinsic to the most powerful radio sources. 

Recall that in the previous subsection, we found that 
stellar population differences between radio-loud galaxies 
and the radio-quiet controls disappears at radio luminosi- 
ties above ~ 10 25 ' 5 W Hz -1 . These highly luminous radio 
galaxies either correspond to a very early phase of the evo- 
lution of the black hole/jet before star formation has been 
suppressed, or they represent a separate population of black 



holes that are undergoing a different accretion mechanism. 
We will investigate these questions in more depth in future 
work. 



4 SUMMARY 

We have cross-matched the SDSS-III BOSS galaxy spectro- 
scopic sample with the NVSS and FIRST surveys, generat- 
ing the largest radio galaxy sample with optical spectroscopy 
at z ~ 0.6 to date. Combining th is sample with th e SPS S 
DR7 radio-loud AGN catalogue of lBest fc Heckmanl (|2012f ). 
we study the evolution of the recent star formation histories 
of the host galaxies of radio- loud AGN, as well as radio-quiet 
"control" galaxies that are matched in redshift, stellar mass, 
and velocity dispersion. 

In previous work (Chen et al 2012), we found that the 
fraction of massive galaxies that have formed a significant 
fraction (> a few percent) of their stars in the past Gyr has 
strongly decreased with redshift from z ~ 0.6 to the present 
day. The goal of this paper is to analyze if this evolution is 
linked to the radio-emitting jets hosted by massive galaxies. 
Our main empirical results may be summarized as follows: 

(i) For galaxies with stellar masses M, > 10 11A M Q , the 
fraction of actively star-forming galaxies is ~ 2 times lower 
among galaxies with radio-emitting jets than in the radio- 
quiet control samples. This is true both at low- and at high- 
redshift. The factor of two difference in the fraction of star- 
forming galaxies is independent of radio luminosity, except 
for radio galaxies with luminosities in excess of 10 25 - 5 W 
Hz -1 , where the difference disappears. 

(ii) The locus of massive galaxies in the D4000— H5A 
plane suggests that their star formation histories are charac- 
terized by bursts rather than low-level continuous star for- 
mation. A smaller fraction of radio-loud galaxies have under- 
gone significant (> a few percent in mass fraction) starbursts 
in the last ~Gyr compared to radio-quiet objects. 

(iii) There is a strong correlation between [O ill] EQW 
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and stellar population age (D4000) both for radio-loud AGN 
and for the control sample of radio-quiet galaxies. The rela- 
tion between [O ill] EQW and D4000 does not change with 
redshift. 

(iv) At fixed radio power, there is strong positive evo- 
lution in [O ill] EQW with redshift. The [O m]/H/3 ratio 
also increases at higher redshift, i.e.; radio galaxies are more 
LINER-like at low redshift and more Seyfert-like at high 
redshift. We interpret these correlations as induced trends, 
which follow from the fact that all high redshift massive 
galaxies have more gas, younger stellar populations and 
larger black hole accretion rates. 

Our results suggest a picture in which massive galaxies 
experience cyclical episodes of gas accretion, star formation 
and black hole growth, followed by production of a radio jet 
that acts to shut down the star formation in the galaxy. The 
behaviour of galaxies with M* > 10 11 Mq is the same at 
2 = 0.6 as it is at z = 0.2, except that higher redshift galax- 
ies experience more black hole growth and star formation 
and produce more luminous radio jets during each accretion 
cycle. 

Let us imagine a scenario in which gas is gradually ac- 
creting onto a massive galaxy. There may be an "accumu- 
lation phase" where the cold gas gradually builds up in a 
disk with little or no star formation. After some critical gas 
density is reached, th e disk may becom e unstable and lose 
angular momentum (|Salim et alj 12012] ). Molecular clouds 
form and star formation is triggered. This corresponds to 
the "burst phase". The radio source is triggered when gas 
finally reaches the central black hole. As the jet expands, it 
pushes away the surrounding cold gas, causing the starburst 
to shut down. The radio lobes expand into the surrounding 
intra-cluster medium, heat the ambient gas and shut down 
cooling onto the central galaxy. The shutdown of cooling re- 
moves the power source of the jet and the galaxy eventually 
returns to its original radio-quiet accumulation phase. 

This picture is consistent with a population of ra- 
dio galaxies with lower cold gas content and older stel- 
lar populations than radio-quiet galaxies. The strong evo- 
lution in the fraction of massive galaxies with significant 
star formation may be explained by increased cooling ef- 
fic iencies in massive hal os at higher redshift (see Figure 2 
of IWhite fc FrenkHl99ll ). If the accumulation and starburst 
phases have roughly equal durations, this would lead to a 2:1 
ratio of starbursts in the radio-quiet and radio-loud popula- 
tions, as observed. Recall, however, that clustering analyses 
demonstrate that radio-loud AGN are pref erentially located 
at the centers of m as sive dark matter halos jBest et al.|[2007l; 
Wake et al]|2008al lbl; iMandelbaum et al.ll2009l ; lDonoso et all 
where gas cooling and accumulation rates are presum- 
ably higher. A fair assessment of the jet/burst duty cycle 
requires one to control for such environment-dependent ef- 
fects. However, in CMASS samples, most galaxies are clas- 
sified as central galaxies, satellites are a very small fraction 
of the total — this is because of the very high stellar mass 
selection (Kovac et al. 2012, in preparation). 

The missing link in our analysis of the accre- 
tion/starburst/jet cycle in massive galaxies is direct observa- 
tions of the gas itself. An analysis of Chandra X-ray Obser- 
vatory archival data of the hot gas in 222 nearby galaxy clus- 
ters revealed that Ha and radio emission from the brightest 



cluster galaxy is much mor e pronounced when the the clus- 
ter's core gas entropy is low (|Cavagnolo et al1l2008h . System- 
atic studies of atomic and molecular gas in massive galaxies 
are still limited to small samples fe.g.lSchawinski et al.ll20"09l ; 
ICatinella et al.l 120101 ; ISaintonge et al.l l201lf ) . Because radio 
galaxies are rare, large wide-field surveys will be required be- 
fore samples are large enough to study how the cold interstel- 
lar medium of massive radio-loud galaxies differs from that 
of their radio-quiet counterparts. In the shorter term, it will 
possible to stack DR7 and BOSS spectra to search for signa- 
tures of ionized gas outflows in luminous radio -loud galax- 
ies a nd in post-burst radio-quiet galaxies (e.g. IChen et al.l 
2010). This will be the subject of future work. 



ACKNOWLEDGEMENTS 

We are very grateful to P. N. Best for kindly providing 
his codes for matching SDSS, NVSS and FIRST. The re- 
search is supported by the National Natural Science Foun- 
dation of China (NSFC) under NSFC-10878010, 10633040, 
11003007 and 11133001, the National Basic Research Pro- 
gram (973 program No. 2007CB815405) and the National 
Science Foundation of the United States Grant No. 0907839. 

Thanks go to the Aspen Center for Physics, which is 
supported by NSF Grant No. 1066293, for hospitality during 
the writing of this paper. 

Funding for SDSS-III has been provided by the Al- 
fred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, and the U.S. Depart- 
ment of Energy. SDSS-III is managed by the Astrophysi- 
cal Research Consortium for the Participating Institutions 
of the SDSS-III Collaboration including the University of 
Arizona, the Brazilian Participation Group, Brookhaven 
National Laboratory, University of Cambridge, University 
of Florida, the French Participation Group, the German 
Participation Group, the Instituto de Astrofisica de Ca- 
narias, the Michigan State/Notre Dame/JINA Participation 
Group, Johns Hopkins University, Lawrence Berkeley Na- 
tional Laboratory, Max Planck Institute for Astrophysics, 
New Mexico State University, New York University, Ohio 
State University, Pennsylvania State University, University 
of Portsmouth, Princeton University, the Spanish Partic- 
ipation Group, University of Tokyo, University of Utah, 
Vanderbilt University, University of Virginia, University of 
Washington, and Yale University. 

The research uses the NVSS and FIRST radio surveys, 
carried out using the NRAO VLA: NRAO is operated by 
Associated Universities Inc., under co-operative agreement 
with the National Science Foundation. 



REFERENCES 

Abazajian K. N., Adelman-McCarthy J. K., Agueros M. A., 
Allam S. S., Allende Prieto C, An D., Anderson K. S. J., 
Anderson S. F., Annis J., Bahcall N. A., et al. 2009, ApJS, 
182, 543 

Aihara H., Allende Prieto C, An D., Anderson S. F., 
Aubourg E., Balbinot E., Beers T. C, Berlind A. A., et 
al. 2011, ApJS, 193, 29 



© 0000 RAS, MNRAS 000, 000-000 



The link between radio AGN activity and star formation 13 



Becker R. H., White R. L., Helfand D. J., 1995, ApJ, 450, 
559 

Best P. N., Heckman T. M., 2012, MNRAS, 421, 1569 
Best P. N., Kaiser C. R., Heckman T. M., Kauffmann G., 

2006, MNRAS, 368, L67 
Best P. N., Kauffmann G., Heckman T. M., Ivezic Z., 

2005a, MNRAS, 362, 9 
Best P. N., Kauffmann G., Heckman T. M., Brinchmann 

J., Chariot S., Ivezic Z., White S. D. M., 2005b, MNRAS, 

362, 25 

Best P. N., von der Linden A., Kauffmann G., Heckman 
T. M., Kaiser C. R., 2007, MNRAS, 379, 894 

Binney J., Tabor G., 1995, MNRAS, 276, 663 

Birzan L., McNamara B. R., Nulsen P. E. J., Carilli C. L., 
Wise M. W., 2008, ApJ, 686, 859 

Birzan L., Rafferty D. A., McNamara B. R., Wise M. W., 
Nulsen P. E. J., 2004, ApJ, 607, 800 

Boehringer H., Voges W., Fabian A. C, Edge A. C, Neu- 
mann D. M., 1993, MNRAS, 264, L25 

Bower R. G., Benson A. J., Malbon R., Helly J. C, Frenk 
C. S, Baugh C. M., Cole S., Lacey C. C, 2006, MNRAS, 
370, 645 

Brinchmann J., Chariot S., White S. D. M., Tremonti C, 
Kauffmann G., Heckman T., Brinkmann J., 2004, MN- 
RAS, 351, 1151 

Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 

Catinella B., Schiminovich D., Kauffmann G., Fabello S., 
Wang J., Hummels C, Lemonias J., Moran S. M., et al. 
2010, MNRAS, 403, 683 

Cavagnolo K. W., Donahue M., Voit G. M., Sun M., 2008, 
ApJ, 683, L107 

Chen Y.-M., Kauffmann G., Tremonti C. A., White S., 
Heckman T. M., Kovac K., Bundy K., Chisholm J., et 
al. 2012, MNRAS, 421, 314 

Chen Y.-M., Tremonti C. A., Heckman T. M., Kauffmann 
G., Weiner B. J., Brinchmann J., Wang J., 2010, AJ, 140, 
445 

Chen Y.-M., Wild V., Kauffmann C, Blaizot J., Davis M., 
Noeske K., Wang J.-M., Willmer C, 2009, MNRAS, 393, 
406 

Churazov E., Briiggen M., Kaiser C. R., Bohringer H., For- 

man W., 2001, ApJ, 554, 261 
Colless M., Dalton G., Maddox S., Sutherland W., Norberg 

P., Cole S., Bland-Hawthorn J., Bridges T., et al. 2001, 

MNRAS, 328, 1039 
Condon J. J., Cotton W. D., Broderick J. J., 2002, AJ, 124, 

675 

Condon J. J., Cotton W. D., Greisen E. W., Yin Q. F., 
Perley R. A., Taylor G. B., Broderick J. J., 1998, AJ, 115, 
1693 

Croton D. J., Springel V., White S. D. M., Dc Lucia G., 
Frenk C. S., Gao L., Jenkins A., Kauffmann G., et al. 
2006, MNRAS, 365, 11 

Donoso E., Best P. N., Kauffmann G., 2009, MNRAS, 392, 
617 

Donoso E., Li C, Kauffmann G., Best P. N., Heckman 

T. M., 2010, MNRAS, 407, 1078 
Dunlop J. S., Peacock J. A., 1990, MNRAS, 247, 19 
Eisenstein D. J., Weinberg D. H., Agol E., Aihara H., Al- 

lende Prieto C, Anderson S. F., Arns J. A., Aubourg E., 

Balbinot E., et al. 2011, AJ, 142, 72 
Fabian A. C, Sanders J. S., Allen S. W., Crawford C. S., 



Iwasawa K., Johnstone R. M., Schmidt R. W., Taylor 

G. B., 2003, MNRAS, 344, L43 
Fabian A. C, Sanders J. S., Ettori S., Taylor G. B., Allen 

S. W., Crawford C. S., Iwasawa K., Johnstone R. M., et 

al. 2000, MNRAS, 318, L65 
Fabian A. C, Sanders J. S., Taylor G. B., Allen S. W., 

2005, MNRAS, 360, L20 

Forman W., Nulsen P., Heinz S., Owen F., Eilek J., 
Vikhlinin A., Markevitch M., Kraft R., et al. 2005, ApJ, 
635, 894 

Fukugita M., Ichikawa T., Gunn J. E., Doi M., Shimasaku 

K., Schneider D. P., 1996, AJ, 111, 1748 
Gunn J. E., Carr M., Rockosi C, Sekiguchi M., Berry K., 

Elms B., de Haas E., Ivezic Z., et al. 1998, AJ, 116, 3040 
Gunn J. E., Siegmund W. A., Mannery E. J., Owen R. E., 

Hull C. L., Leger R. F., Carey L. N., Knapp G. R., et al. 

2006, AJ, 131, 2332 

Heckman T. M., Kauffmann G., Brinchmann J., Chariot 

S., Tremonti C, White S. D. M., 2004, ApJ, 613, 109 
Hogg D. W., Finkbeiner D. P., Schlegel D. J., Gunn J. E., 

2001, AJ, 122, 2129 

Johnston H. M., Sadler E. M., Cannon R., Croom S. M., 

Ross N. P., Schneider D. P., 2008, MNRAS, 384, 692 
Karim A., Schinnerer E., Martmez-Sansigre A., Sargent 

M. T., van der Wei A., Rix H.-W., Ilbert O., Smolcic V., 

et al. 2011, ApJ, 730, 61 
Kauffmann G., Heckman T. M., Best P. N., 2008, MNRAS, 

384, 953 

Kauffmann G., Heckman T. M., White S. D. M., Chariot 

S., Tremonti C, Brinchmann J., Bruzual G., Peng E. W., 

et al. 2003a, MNRAS, 341, 33 
Kauffmann G., Heckman T. M., White S. D. M., Chariot 

S., Tremonti C, Peng E. W., Seibert M., Brinkmann J., 

et al. 2003b, MNRAS, 341, 54 
Kauffmann G., Heckman T. M., Tremonti C, Brinchmann 

J., Chariot S., White S. D. M., Ridgway S. E., Brinkmann 

J., et al. 2003c, MNRAS, 346, 1055 
Mandelbaum R., Li C, Kauffmann C, White S. D. M., 

2009, MNRAS, 393, 377 
Mauch T., Sadler E. M., 2007, MNRAS, 375, 931 
Padmanabhan N., Schlegel D. J., Finkbeiner D. P., Baren- 

tine J. C, Blanton M. R., Brewington H. J., Gunn J. E., 

Harvanek M., et al. 2008, ApJ, 674, 1217 
Pier J. R., Munn J. A., Hindsley R. B., Hennessy G. S., 

Kent S. M., Lupton R. H., Ivezic Z., 2003, AJ, 125, 1559 
Rawlings S., Jarvis M. J., 2004, MNRAS, 355, L9 
Rola C. S., Terlevich E., Terlevich R. J., 1997, MNRAS, 

289, 419 

Sadler E. M., Jackson C. A., Cannon R. D., Mclntyre V. J., 
Murphy T., Bland-Hawthorn J., Bridges T., Cole S., et al. 

2002, MNRAS, 329, 227 

Saintonge A., Kauffmann G., Kramer C, Tacconi L. J., 

Buchbender C, Catinella B., Fabello S., Gracia-Carpio 

J., et al. 2011, MNRAS, 415, 32 
Salim S., Fang J. J., Rich R. M., Faber S. M., Thilker D. A., 

2012, ArXiv e-prints 
Schawinski K., Khochfar S., Kaviraj S., Yi S. K., Boselli 

A., Barlow T., Conrow T., Forster K., et al. 2006, Nature, 

442, 888 

Schawinski K., Lintott C. J., Thomas D., Kaviraj S., Viti 
S., Silk J., Maraston C, Sarzi M., et al. 2009, ApJ, 690, 
1672 



© 0000 RAS, MNRAS 000, 000-000 



14 Chen et al. 



Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Strauss M. A., Weinberg D. H., Lupton R. H., Narayanan 

V. K., Annis J., Bernardi M., Blanton M., Buries S., et 

al. 2002, AJ, 124, 1810 
Tremonti C. A., Heckman T. M., Kauffmann G., Brinch- 

mann J., Chariot S., White S. D. M., Seibert M., Peng 

E. W., et al. 2004, ApJ, 613, 898 
Wake D. A., Sheth R. K., Nichol R. C, Baugh C. M., 

Bland-Hawthorn J., Colless M., Couch W. J., Croom 

S. M., et al. 2008a, MNRAS, 387, 1045 
Wake D. A., Croom S. M., Sadler E. M., Johnston H. M., 

2008b, MNRAS, 391, 1674 
White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 
York D. C, Adelman J., Anderson Jr. J. E., Anderson 

S. F., Annis J., Bahcall N. A., Bakken J. A., Barkhouser 

R., et al. 2000, AJ, 120, 1579 
Zheng X. Z., Bell E. F., Papovich C, Wolf C, Meisen- 

heimer K., Rix H.-W., Rieke G. H., Somerville R., 2007, 

ApJ, 661, L41 



© 0000 RAS, MNRAS 000, 000-000 



